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ABSTRACT 


Objectives: Since soil carbon (C) is a principal source of energy for the nutrient-recycling activities of heterotrophic soil organ¬ 
isms, the maintenance of belowground C stocks is vital for sustaining forest productivity. 

Methods: The present study was conducted in temperate region of the Garhwal Himalaya during the year 2014-2016 to assess 
the belowground carbon flux in the forest ecosystem. In the study, we studied the component wise belowground carbon flux in 
trees, soil organic carbon (SOC) and litter carbon of six different forest types for measuring total belowground carbon allocation 
(TBCA). 

Results and Conclusion: The total belowground carbon of live trees varies 20.02 to 60.58 MgC/ha, whereas stock root carbon 
(14.01-38.27 MgC/ha), lateral roots carbon (5.24-17.57 MgC/ha) and fine root carbon (0.67-12.2 MgC/ha) in selected forest 
types were recorded. The maximum SOC was exhibited by Abies pindrow forest (110.83± 5.04 MgC/ha), followed by Pinus rox- 
burghii forest (108.22±13.03MgC/ha), Quercus floribunda forest (97.37±7.64 MgC/ha), whereas minimum SOC was recorded 
in Cedrus deodara forest (56.94±5.13 MgC/ha). The maximum value of litter carbon was recorded for Abies pindrow forest 
(2.94±1.02 MgC/ha), followed by Quercus semecarpifolia forest (2.22±0.33 MgC/ha), Quercus floribunda forest (2.06±0.28 MgC/ 
ha), Cedrus deodara forest (1.86±0.26 MgC/ha), Quercus leucotrichophora forest (1.44±0.27 MgC/ha), Pinus roxburghii forest 
(0.84±0.10 MgC/ha). Forest ecosystem in Garhwal Himalaya appears to be the most conducive soil-climatic environment for 
higher accumulation of SOC, thus helping in maintaining the soil quality. The study showed that belowground carbon stocks in 
Abies pindrow forests seems has maximum in carbon assimilatory capacity, whereas Cedrus deodara forest has minimum BGC 
stocks. There is a huge potential to increase SOC potential through the soil conservation and hence should be implemented. 
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INTRODUCTION 

Carbon sequestration through forestry has immense potential 
and plays a significant role in solving critical global envi¬ 
ronmental problems such as atmospheric accumulation of 
Green House gases (GHG) and climate change. Estimates 
of existing C stock pools, stored in various forest types can 
be helpful in making decisions about C management. For¬ 
est biomass constitutes the largest terrestrial carbon sink 
for CO, removal from the atmosphere through the process 
of photosynthesis. Forests absorb CO, from the surrounding 
atmosphere and store carbon in their different components 


like wood, leaves, litter, roots and soil allocating as “carbon 
sinks” (Haripriya, 2003). The major carbon pools in forest 
ecosystem are plant biomass (above and below ground), 
coarse woody debris, litter and soil (Sharma et ah, 2016a). In 
most forested ecosystems, the majority of the carbon is stored 
below ground as roots and decaying biomass or as dead or¬ 
ganic matter (DOM). Forest floor contribute large amounts 
of organic material to the soil in the form of different com¬ 
ponents such as leaves, twigs, branches, reproductive parts, 
fruits, where their decomposition releases different nutrients 
into the soil (Tandel et ah, 2009). However, belowground 
biomass (stock, lateral and fine roots components of trees) 
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also plays an important role in carbon sequestration. Total 
Belowground Carbon (TBGC) is a large fraction of gross pri¬ 
mary production (more than 30%) (Ryan et al., 1994; Gower 
et al., 1996). In some ecosystem, amount of TBGC can sur¬ 
pass the values of aboveground net primary production (Law 
et al., 1999) and provides the principal source of detritus C 
to mineral soil. Despite the significant role of TBGC in the 
C budget of terrestrial ecosystems, controls on TBGC are 
poorly understood. There is a large dissimilarities in the rate 
and the length of time that forest floor carbon may seques¬ 
tered in soil that are related to the vegetation productivity, 
biological and physical conditions in the soil, the past history 
of soil organic carbon inputs and various disturbances in¬ 
cluding physical and anthropogenic (Post and Kwon, 2000). 

Forest ecosystem in Garhwal Himalaya appears to be the 
most conducive soil-climatic environment for higher accu¬ 
mulation of SOC, thus helping in maintaining the soil qual¬ 
ity. It is essential to examine the changes in carbon fluxes 
derived from land-use change patterns to obtain basic infor¬ 
mation on the carbon contents associated with the various 
stocks of the natural forests. In earlier studies from Garh¬ 
wal Himalaya, various researchers (Chaturvedi et al., 1982; 
Rawat and Singh, 1988; Adikari et al., 1995) from Kumoun 
Himalaya and (Sharma et al., 2010; Gairola et al., 2011; 
Sharma et al., 2016c) from Garhwal Himalaya attempted to 
predict biomass and carbon stocks. The present paper aimed 
at the assessment of the forest diversity, componentised con¬ 
tribution of belowground carbon accumulation of tree spe¬ 
cies, litter fall and soil organic carbon (SOC) stocks of the 
different forests of Garhwal Himalaya under temperate con¬ 
ditions. 


METHOD AND METHODOLOGY 

Study area 

The Uttarakhand state is located in North-West part of the 
country. Its geographical area is 53,483 km 2 which consti¬ 
tutes 1.63% of total area of the country (FSI, 2013). The pre¬ 
sent study was conducted in temperate region of the Garhwal 
Himalaya during the year 2014-2016 to assess the below¬ 
ground carbon flux in the forest ecosystem. The study area is 
located in Uttarkashi and Tehri district of Garhwal Himalaya 
along the catchment area of Bhagirathi River (tributary of 
river Ganga). In this study, general survey of the area was 
carried out at selected sites of Tehri and Uttarkashi district of 
Garhwal Himalaya (Figurel). After the reconnaissance sur¬ 
vey, we selected following six dominant forest types were 
named according to the classification given by Champion 
and Seth (1968) viz, 

(i) Finns wallichiana forest (PF) - Pine Forest (9/Clb) 

(ii) Quercus leucotrichophora forest (QLF) - Ban Oak 
forest (12/Cla) 


(iii) Quercus floribunda forest (QFF) - Mom Oak Forest 
(12/Clb) 

(iv) Quercus semecarpifolia forest (QSF) - West Himala¬ 
yan Upper Oak Forest (12/C2b) 

(v) Cedrus deodara forest (CF) - Dry Deodar Forest (13/ 
C2b) 

(vi) Abies spectabilis forest (AF) - West Himalayan Sub- 
alpine Fir Forest (14/Cla) 

Methodology 

To analyse the forest vegetation on different ridge tops, 10 
sample plots of 0.1 ha each were laid out in 6 selected ridge 
tops of each forest types (10 plots x 06 forest types = total 
60 sample plots). All individuals > 10 cm diameter at breast 
height (DBH = 1.37 m from ground) were considered as tree 
in each sample plot. The DBH and height of all the trees fall¬ 
ing within the sample plot were measured by tree Calliper 
and Ravi multimeter respectively. The tree height on differ¬ 
ent slope positions was measured following MacDicken et 
al. (1991). The slope correction was employed for the sample 
plots which are located on a slope > 10%, so that the adjust¬ 
ment can be made to the plot area at the time of analysis. The 
slope angle was measured by clinometer. 

The biomass of the tree species were calculated by regression 
equations. The tree components (bole, bole bark, branches, 
twigs, leaves, stump roots, lateral roots, fine roots) were cal¬ 
culated by various equations developed by Rawat and Singh 
(1988), Garkoti and Singh (1992) and Adhikari et al. (1995). 
Biomass equations by tree components were developed to 
relate oven dry weight to tree cbh. The fonn of the allometric 
function of the equation was: 

LnY=a + ALnX 

Where, Y is weight of tree component in kilograms, and X is 
tree circumference in centimetres measured at breast height 
(CBH), a and b is the intercept and slope of the particular tree 
species respectively. The total Carbon density (TCD) was es¬ 
timated by the following formula (Sharma et al., 2010): 

Carbon (C Mg/ha) = Biomass (Mg/ha) xCarbon % 

Forest litter was collected by using I x | m randomly placed 
quadrat at three places in 0. Ilia sample plot. Litter was col¬ 
lected three times in a year consisted of basically fresh and 
partially decomposed leaves, bark and reproductive parts. 
Fresh weight was determined in the field. The collected litter 
was brought to the laboratory and oven dried at 80°C up to 
the constant weights. The biomass hence obtained from the 
detritus was then multiplied by appropriate carbon fraction 
according to IPCC (2006) and extrapolated for a hectare. 

For organic C determination, the soil samples were sieved 
through a 2 mm sieve and then thoroughly mixed. Modified 
Walkley and Black’s rapid titration method (Walkley, 1947) 
was used to estimate the SOC content in the collected soil 
samples (Mehta et al., 2014). The contents of organic carbon 
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in soil estimated in percentage were then converted to tonnes 
per hectare using bulk density, depth of soil and area. 

soil mass in layer X SOC concentration (%) 
SOC (MgC/ha) =-—-— 


RESULTS 

The mean values of density, species richness and below¬ 
ground carbon allocation in different forest types are shown 
in table 1. A total of 27 tree species were occurred in all six 
studied forest types, out of which highest 12 species were 
recorded in QFF and least 8 species were found in CF. The 
highest stem density was seen AF (778 ± 89.97 trees/ha) fol¬ 
lowed by QFF and the lowest density seen in pine forest (560 
± 69.93 trees/ha). The different studied forests showed their 
distribution above 1500m asl. upto 3500 m asl. with QFF 
have showed narrow distributional range whereas, QLF and 
QSF distributed over broad elevational range. 

The total live tree carbon stocks in these studied forest types 
were seen in the range of 163.47 ± 14.08 Mg C ha' 1 to 320.3 
± 29.7 Mg C ha 1 , Out of which the variation in TBCA seen in 
the range of 20.02 ± 2.47 to 60.58 ± 5.93 Mg C ha 1 . AF had 
the highest value of carbon stock followed by PF, CF had the 
lowest TBCA in case of conifer dominated forests. In case of 
broad-leaved forests, QSF forests has highest carbon stocks 
value, whereas QFF and QLF had showed carbon stocks in 
the range of 49.96 ± 5.8 Mg C ha' 1 and 47.07 ± 5.58 Mg C 
ha -1 respectively. In our study, the component wise distribu¬ 
tion of carbon in burial parts of trees shows that stock roots 
(14.01-38.27 MgC/ha) stores more carbon than other roots 
followed by lateral roots carbon (5.24-17.57 MgC/ha) and 
fine root carbon (0.67-12.2 MgC/ha) in selected forest types 
were recorded. In different forest types, QLF shows highest 
value of Stock roots Carbon, AF shows highest value of Lat¬ 
eral roots carbon and QFF shows highest carbon storage in 
lateral roots system. However, mean value of aboveground 
biomass and carbon stocks are seen in highest in AF (259.71 
± 23.87 MgC/ha) followed by QFF (215.23 ± 20.58 MgC/ha) 
and least value in PF (139.11 ± 11.44 MgC/ha). 

The soil organic carbon (SOC) in different forest types 
are shown in Table 1. The maximum value of litter carbon 
was recorded for Abies pindrow forest (2.94±1.02 MgC/ 
ha), followed by Quercus semecarpifolia forest (2.22±0.33 
MgC/ha), Quercus floribunda forest (2.06±0.28 MgC/ha), 
Cedrus deodara forest (1.86±0.26 MgC/ha), Quercus leu- 
cotrichophora forest (1.44±0.27 MgC/ha), Pinus roxburghii 
forest (0.84±0.10 MgC/ha). The maximum SOC was ex¬ 
hibited by Abies pindrow forest (110.83± 5.04 MgC/ha), 
followed by Pinus roxburghii forest (108.22±13 MgC/ha), 
Quercus floribunda forest (97.37±7.64 MgC/ha), where¬ 
as minimum SOC was recorded in Cedrus deodara forest 
(56.94±5.13 MgC/ha). The correlation between different 


carbon components and diversity indices are shown in ta¬ 
ble 2. Total basal cover and species richness show negative 
correlation, however density shows positive correlation with 
AGC and TCD. The patterns of variations in different forest 
types across different ecological varibles are presented by 
mean of PCA diagram (Figure 2). The distance between eco¬ 
logical parameters (blue dots) approximates the dissimilarity 
of distribution of relative abundance of those species across 
these forests. The distance between forests (red dots) show 
which ecological parameter effects the forest composition 
within the forest and between the forests. 


DISCUSSION 

Along with the climatic variations, lower elevational QLF 
had comparatively higher number of species than lower 
number of species was found at higher elevational AF, which 
implies the climatic adaptation of plant species. There is also 
seen exponential decline in biomass and carbon increment 
with increasing DBH as observed in our study is primarily 
related to the age of trees (Figure 3). Elevation and tempera¬ 
ture, rainfall are main physical factors for forest structure 
and carbon allocation in different in Garhwal Himalaya. For¬ 
ests in studied area were mature with higher girth values as 
they were undisturbed. According to Saxena et al. (1979), 
trees with higher girths indicate the best representation of 
a species in the particular forest in specific environmental 
conditions whereas, lower girths either indicate the chance 
occurrence of the species in that area or show presence of the 
biotic disturbance in the past. Age-related declines in carbon 
stocks are widely documented (e.g., Pregitzer and Euskirch- 
en, 2004; Bradford and Kastendick, 2010; D’Amato et al., 
2011 ). 

Tree biomass and carbon stocks in forest ecosystems vary 
with forest type, species composition, stand age, size class 
of trees, site conditions, rainfall pattern and altitude (Sharma 
et al., 2016a; 2016b; 2016c; Gairola et al., 2011; Zhao et al., 
2014). The values of belowground biomass density in the 
present study ranged between of 20.02 ± 2.47 to 60.58 ± 5.93 
Mg C ha'. Negi et al. (2003) observed that the conifers have 
maximum C stored followed by mixed and broad leaved for¬ 
ests. The BGCA values obtained were lowered than Sheikh 
et al 2012 but followed the values reported by Sharma et al. 
(2010). These values are higher than earlier reported values 
from Garhwal Himalaya, comparative values of TBD and 
TCD from Uttarakhand and other parts of India. 

Type stratum AF in general exhibited higher values for high¬ 
er SOC density (110.83± 5.04 MgC/ha) followed by in PF 
(108.22±13.03 MgC/ha) and lowest value of SOC were cal¬ 
culated in QSF (74.65±6.10 MgC/ha) which may be due to 
the presence of mature girth classes compared to other type 
strata and existence of coniferous leaves, which generally 

n 
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decompose slower than leaves of broadleaf species (Mendo¬ 
za-Vega et al., 2003). Lower temperatures may be also one 
of the reasons for highest carbon stocks at higher altitudes 
Schlesinger (1997). However, Liebens and VanMolle (2003) 
also reported the soils of coniferous forest might store more 
SOC per hectare than broadleaf forest soils. Litter fall carbon 
also shower highest value in AF (2.94 ± 1.02 MgC/ha) and 
lowest in PF (0.84 ±0.1 MgC/ha). Hobbie et al. (2006) also 
found that litter turnover rate was positively correlated with 
mean annual soil temperature. 

The study showed that belowground carbon stocks in Abies 
pindrow forests seems has maximum in carbon assimilatory 
capacity, whereas Cedrus deodam forest has minimum BGC 
stocks. These Forest floor thus play a key role in the global 
carbon budget in Garhwal Himalaya and can have large im¬ 
pact on carbon release under a climate change scenario (Lai, 
2002). These forests have huge potential to increase SOC 
potential through the soil conservation and hence should be 
implemented. 

As the climate change issues became prominent on political 
and corporate agenda, it is duty of people of India and other 
countries to start recognizing their responsibility towards 
taking action against global warming. The prevention of de¬ 
forestation and promotion of afforestation have often been 
cited as strategies to slow down global warming and climate 
change (Bala et al., 2007) and help in increase to sinks more 
carbon from their present potential. Agro forestry practices 
are the agents that can enhance the ability of forest to sink 
more carbon (Ahmed et al., 2016). Alternate production of 
energy like Hydro power and solar energy can significantly 
reduce the pressure on forest in terms of fuel wood removal. 

CONCLUSION 

The present study showed the AF, CF and oak forests and 
stored more biomass and carbon stocks than other forests; 
hence these forests have higher potential C sinks to tone 
down the global warming consequences as well as meet our 
future energy demands. Data obtained through regular as¬ 
sessments of these carbon stocks will yield knowledge on the 
impacts of particular conservation and management regimes 
on the forest resources that can assist the enhanced carbon 
management. Optimized forest management with regard to 
conservation should implement to secure a high productiv¬ 
ity of the forest and avoid disturbances as much as possible. 
Sustainable use of the forest resources will enable us to con¬ 
serve them for future generation. 
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Table i: Component wise allocation of carbon stocks in different forest types 


Parameters 

Abies 

spectabilis 

Cedrus 

deodara 

Forest Types 

Pinus Q.floribunda 

wallichiana 

Q- 

leucotrichophora 

Q. semecarpifolia 

Species richness 

io 

8 

11 

12 

IO 

IO 

Density 

778 ± 89.97 

592 ± 94.41 

560 ± 69.93 

662 ± 18.81 

616 ± 40.69 

574 ± 45-56 

TBC 

20.73 

21.18 

14.69 

14.60 

18.50 

17.88 

AGC (Mg/ha) 

259.71 ± 23.87 

148.4 ± 27.12 

i39.ll ± 11.44 

215.23 ± 20.58 

162.75 ± 21.32 

174.22 ± 18.48 

SRC (Mg/ha) 

35-84 ± 3-7 

14.01 ± 1.62 

17.99 ± t-86 

25.I9 ± 1.12 

38.27 ± 4.64 

31.24 ± 3.47 

LRC (Mg/ha) 

17.57 ± i-75 

5.24 ± 0.72 

5.7 ± 0.68 

12.56 ± 0.59 

7.54 ± 0.93 

13.65 ± 1.42 

FRC (Mg/ha) 

7.18 ± 0.61 

O.77 ± 0.22 

0.67 ± 0.17 

12.2 ± 7.16 

1.27 ± 0.5 

10.79 ± 5-45 

BGC (Mg/ha) 

60.58 ± 5.93 

20.02 ± 2.47 

24.36 ± 2.65 

49.96 ± 5.8 

47.07 ± 5.58 

55.67 ± 3.88 

TCD (Mg/ha) 

320.3 ± 29.7 

168.42 ± 27.19 

163.47 ± 14 08 

265.18 ± 18.4 

209.82 ± 26.82 

229.89 ± 19.52 

Litter (Mg/ha) 

2.94 ± 1.02 

1.86 ± 0.26 

0.84 ± 0.1 

2.06 ± 0.28 

1.44 ± 0.27 

2.22 ± O.33 

SOC (Mg/ha) 

110.83 ± 5-°4 

56.94 ± 5.13 

108.22 ± 13.03 

97.37 ± 7.64 

80.70 ± 8.32 

74.65 ± 6.10 
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Table 2: Correlation table between ecological variables studied in different forest types 



Species richness 

Density 

TBC 

AGC 

SRC 

LRC 

FRC 

BGC 

TCD 

Litter SOC 

Species richness 

1 










Density 

0.130 

1 









TBC 

-.869* 

0-354 

1 








AGC 

0.294 

■ 947 ** 

0.194 

1 







SRC 

0.204 

0.502 

0.173 

0-535 

1 






LRC 

0.302 

0-754 

0.135 

.906* 

0.620 

1 





FRC 

0.515 

0.349 

-0.267 

0.627 

0.320 

0.785 

1 




BGC 

0.370 

0.625 

0.056 

0.776 

.865 

.903* 

0-733 

1 



TCD 

0.329 

.901* 

0.164 

.984** 

0.653 

.948” 

0.686 

.876* 

1 


Litter 

-0.138 

0-775 

0-533 

.851* 

o -435 

.865* 

0.625 

0.705 

.850* 

1 

SOC 

0.722 

0.505 

-0.460 

0.519 

0.260 

o -433 

0.199 

0.342 

0-494 

0.055 1 


*. Correlation is significant at the 0.05 level ( 2 -tailed). 
**. Correlation is significant at the 0.01 level ( 2 -tailed). 



Figure 1: Map of the Study Area ( not up to scaled). 
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Figure 3: Forest structure in term of DBH distribution. 


Figure 2: PCA diagram between ecological variables. 
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